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ABSTRACT 

The magnetospheric locations of pulsar radio emission region are not well known. The 
actual form of the so-called radius-to-frequency mapping should be reflected in the 
aberration-retardation (A/R) effects that shift and/or delay the photons depending 
on the emission height in the magnetosphere. Recent studies suggest that in a handful 
of pulsars the A/R effect can be discerned w.r.t the peak of the central core emission 
region. To verify these effects in an ensemble of pulsars we launched a project analysing 
multi-frequency total intensity pulsar profiles obtained from the new observations 
from the Giant Meterwave Radio Telescope (GMRT), Arecibo Observatory (AO) and 
archival European Pulsar Network (EPN) data. For all these profiles we measure 
the shift of the outer cone components with respect to the core component which is 
necessary for establishing the A/R effect. Within our sample of 23 pulsars 7 show the 
A/R effects, 12 of them (doubtful cases) show a tendency towards this effect, while 
the remaining 4 are obvious counter examples. The counter-examples and doubtful 
cases may arise from uncertainties in determination of the location of the meridional 
plane and/or the core emission component. It hence appears that the A/R effects are 
likely to operate in most pulsars from our sample. We conclude that in cases where 
those effects are present the core emission has to originate below the conal emission 
region. 

Key words: pulsars: general - stars: neutron - radiation mechanisms: non-thermal 
- methods: data analysis 

1 INTRODUCTION 

Aberration and retardation effects (A/R hereafter) can be observed in a pulsar profile as a shift of the position of conal 
components with respect to the core component towards earlier longitudinal phases (see for e.g. Malov & Suleimanova (1998), 
Gangadhara & Gupta 2001, G&Ga hereafter). Such effects should occur if different components are emitted at different 
distances from the pulsar surface (emission radii), as well as from the pulsar spin axis. Aberration is caused by bending of 
radiation beam due to the polar cap rotation, while retardation is based on a path difference for radiation from different 
conal emission regions to reach an observer. If we assume that the emission from the core component arises relatively close to 
the star surface, then it should not be strongly affected by either of the two above mentioned effects. This will be our initial 
assumption. 

To determine A/R shifts the pulsar profile has to meet certain requirements. It has to be a high resolution profile with high 
signal to noise (S/N) ratio. The core and the conal components have to be clearly identified within the profile. Multifrequency 
data is recommended, so one can follow the profile evolution throughout all frequencies, which can help to identify different 
emission components. 

When values of A/R shifts are determined, then the heights of the radiation emission region (emission altitudes hereafter) 
can be calculated (see G&Ga and Dyks et. al 2004). It is believed that at different frequencies the emission arises at different 
heights above the pulsar surface (Kijak & Gil 1998, 2003 and Mitra & Rankin 2002). The results of this analysis can be used 
to verify the existence of a radius-to-frequency mapping. All observational limits for emission altitude hence can be crucial 
for understanding the physical mechanism of generation of pulsar coherent radio emission. 
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The relativistic beaming model initially proposed by Blaskiewicz, Cordes & Wasserman (1991, BCW hereafter) clearly 
demonstrated that aberration and retardation effects play an important role in pulsars. This study was primarily based on 
evidence which followed from the effects of A/R as seen in the inflexion point of the pulsar's polarisation position angle 
(hereafter PA) traverse, which lags the midpoint of the total intensity profile centre. A similar effect of A/R was reported by 
G&Ga and Gupta & Gangadhara (2003, G&Gb hereafter) in a handful of pulsars where the core emission was seen to lag 
behind the profile centre. 

In this paper we have undertaken a careful study to establish the A/R effect observed by G&Ga and G&Gb for a large 
sample of pulsars observed at multiple frequencies. Most of the data are new observations from the Giant Meterwave Radio 
Telescope (GMRT hereafter) and the Arecibo Observatory (AO hereafter). We have also used some archival data from the 
European Pulsar Network (EPN hereafter) archiv^. In section ([5)) we discuss various methods used to find emission heights 
in pulsars, in section Q we discuss various factors affecting A/R measurements in pulsars and section Q deals with the 
observation and data analysis methods used in this paper. 

As a result of our analysis presented in section ((5} we found that out of 23 pulsars in our sample 7 clearly show the 
A/R effect, 12 show a clear tendency towards this effect, while the remaining 4 are counter examples. However, as argued in 
section ([3]), all problematic cases (pulsar profiles at all or some frequencies not showing the A/R effect) can be attributed to 
a number of effects like incorrect identification of the core component or missing conal emission. We can conclude that A/R 
effects are seen to operate in pulsars, which we discuss in section ([6)). 



2 GEOMETRY AND EMISSION ALTITUDES 

Radio emission heights in pulsars are primarily obtained by two methods; the geometrical method and heights estimation 
based on A/R effects. Here we briefly mention the essential ingredients of the methods used, and a detailed discussion of the 
various methods used can be found in Mitra & Li (2004). 

Radio emission geometry is determined by several parameters: a — an inclination angle of the magnetic dipole with 
respect to the rotation axis, /3 — the minimum angle between an observer's line of sight and magnetic a:xis (impact angle), p 
— an opening angle of the radio emission beam, r — a radial distance of the radio emission region measured from the centre 
of the neutron star (emission altitude). The opening angle p of the pulsar beam corresponding to the pulse width W is given 
by: 



p = 2 arcsin ^ y sin (a + /3) • sin a ■ sin^ + sin^ ^ J ' 

where a, f3, W and p are measured in degrees (Gil et al. 1984). The opening angle p is the angle between the pulsar magnetic 
axis and the tangent to magnetic field lines at points where the emission corresponding to the apparent pulse width W 
originates. For dipolar field lines: 

p = 1.24°s(r/10''cm)''''p-''''^ , (2) 

(Gil & Kijak 1993), where s — d/rp is a mapping parameter which describes the locus of corresponding field lines on the polar 
cap (s = at the pole and s = 1 at the edge of the polar cap), d is the distance of the given magnetic open field line from the 
dipolar magnetic axis (in cm), Vp is the polar cap radius (in cm) and P is the pulsar period in seconds. The radio emission 
altitude can be obtained using eqn. ([2}: 

'^.-^IOP(^)' km. (3) 

In this equation parameter s = 1 is used which corresponds to the bundle of last open magnetic field lines. Kijak and Gil 
(1997) also derived a semi-empirical formula for emission height which was slightly modified by Kijak & Gil (2003) by using 
larger number of pulsars and broadening the frequency coverage in their analysis. They estimated the emission heights for 
a number of pulsars at several frequencies using the geometrical method (eqs. [T] and ^ based on the above assumption that 
the low intensity pulsar radiation (corresponding to the profile wings) is emitted tangentially to the bundle of the last open 
dipolar field lines. In the method, they used low-intensity pulse width data at about 0.01 % of the maximum intensity (Kijak 
& Gil 1997 and Kijak & Gil 2003). The most recent formula has a form 

//inn _L on^ -0.26±0.09 d0.07±0.03 r)0.30±0.05 , „ 

r-KG = (400 ± 80)!/ghz ^-i5 P km, (4) 

where vghz is observing frequency in GHz and P_i5 = P/10~^^ is the period derivative. The formula tkg sets lower and 
upper boundaries for emission height values in the statistical sense. 



^ |http://www.mpifr-bonn.mpg.de/div/puIsar/data/| 



On the aberration/retardation effects in radio pulsars 3 



For calculating the emission heights one can take a different approach. The A/R effects have been observed in the total 
intensity profiles as discussed by Malov & Suleimanova (1998), G&Ga, G&Gb, Dyks et al (2004). In turn, they have been 
used to calculate radio emission heights. Observationally the A/R effects are obtained by measuring the positions of the conal 
emission components and identifying the central core emission component. The shape of the observed profile depends on the 
structure of radio beam, which is thought to be composed of a central core emission surrounded by nested cones (Rankin 
1983, Mitra & Deshpande 1999). The line of sight passing through the central region of this beam will hence cut the cone 
emission at the edges including the central core region. Following this, pulsar profiles can be classified due to the number 
and type of visible components (Rankin 1983) . If we assume that in a certain profile the core component and both the conal 
components are present, we can define the phase difference: 

A^=^^, (5) 

where cj>i and (j>t are phase positions (in degrees) of the peaks of the leading and trailing components with respect to the core 
phase {4>c = 0°). Following the refinement of G&Ga by Dyks et al. (2004) one can estimate the radio emission altitud^: 

'^^^"-—360^^^^' 

where -Rlc = cP/2-k (where c is speed of light) is the light cylinder radius. The A/R method has an edge over the geometrical 
method since it only depends weakly on the geometrical angles of the star which are difficult to find. Generally the emission 
altitudes found using both these methods are of the order of few hundred kilometers above the surface of the neutron star. In 
this paper we use equation © for estimating emission altitudes following the method outlined by G&Ga. The various factors 
affecting usage of this method are discussed in section ([3} . 



3 FACTORS AFFECTING A/R MEASUREMENTS IN PULSARS 
3.1 Location of the meridional point 

The crucial aspect for estimating radio emission heights using the A/R technique is to identify the location of the meridional 
plane i.e. the plane containing the magnetic axis and the pulsar rotation axis. 

Generally polarisation observations have been used to identify the meridional plane as per the rotating vector model 
(RVM) proposed by Radhakrishnan & Cooke (1969). According to the RVM the variation of the PA of the linear polarisation 
across the pulse reflects the underlying structure of the dipolar magnetic field if one is located on the pulsar frame. In this 
case the steepest gradient or the inflexion point of the PA traverse is the point contained in the meridional plane, and should 
be coincident with the midpoint of the pulse profile. However in the observers frame, as was pointed out by BCW, Dyks 
et al (2004), the A/R effects causes the PA inflexion point to appear delayed w.r.t. the midpoint of the pulse profile. The 
magnitude of the delay is proportional to the height of emission from the neutron star surface and inversely proportional to 
the pulsar period. Hence in the observer's plane the inflexion point of the PA traverse is no longer the point containing the 
meridional plane, rather the meridional plane lies somewhere in between the proflle midpoint and the inflexion point of the 
PA traverse. 

We use the method employed by G&Ga,b to flnd emission heights due to A/R effects using eq. (6). If the mid-point of 
the outer conal peaks leads the core peak, it would imply that the conal emission is higher than the core emission and the shift 
can be used to find emission heights w.r.t. the core. This argument however assumes that the meridional plane passes through 
the peak of the core emission. While this assumption is not rigorously verified, there are few pulsars which have been used 
to justify this argument. One such pulsar is PSR B1857— 26 (see Fig. [l| where the peak of the core, the infiexion point and 
the sign changing circular, are all close to each other (see Mitra & Rankin 2008 for a detailed discussion). If one assumes that 
the core emission is curvature radiation which fills the full polar cap as was suggested by Rankin (1990), then the meridional 
plane would lie in the plane containing the inflexion point, peak of the core and sign changing circular. Although a critical 
look at PSR B1857— 26 reveals that the peak of the core is delayed by a degree or so w.r.t. the inflexion point (see Mitra & 
Rankin 2008). Mitra & Li (2004) have shown that for a few other pulsars the peak of the core lags the inflexion point by 
about a degree. The shifts expected due to A/R effects for slow pulsars are also often small (e.g. for emission altitude of 300 
km and a 0.7 sec pulsar, the expected shift is of the order of a degree) and is comparable to core shifts seen w.r.t. the inflexion 
point. Hence, it is important to note that to demonstrate the A/R effects as suggested by G&Ga and G&Gb, the meridional 



^ Recently in a detailed study Gangadhara (2005) has obtained a revised expression (their eqn. 34) for the relativistic phase shift which 
depends on the angles a and (3, and the corrections in obtaining emission heights turns out to be of the order of 2% to 10% for normal 
pulsars. We have however not used this formula since all the pulsars in our sample are normal pulsars and, further, the error in estimating 
emission heights are often more than 10%. 
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plane should pass through the peak of the core emission (which is currently not verified) and the conal emission should be 
rising higher than the core emission. 

3.2 Identifying the Core Emission 

Identification of the core emission is not trivial in pulsars. Under the core-cone model of pulsar beaming the core is centrally 
located and surrounded by nested cones. This implies that the core is located in the centre of a pulse profile, which is very 
close to the infiexion point of the PA. Since the impact parameter needs to be small to cut through the core emission region 
in a pulsar, the corresponding PA traverse as per the RVM should show a significant rotation across the pulse. Core emission 
is often associated with a sign changing circular polarisation and has a steeper spectral index than the conal emission. 

There are situations where several core like emission characteristics are seen in a pulse component. For example in Fig 
PSR B1831— 04 has a bright central component which is identified as core emission by Rankin (1993a, 1993b). This component 
has a steep spectral index and a sign changing circular as seen in Gould & Lyne's (1998, GL98 hereafter) multi-frequency 
polarisation profiles. On the other hand the PA traverse is fiat, showing no significant rotation (except 90° jumps) and the 
profile is extremely wide (~ 150°). Such emission properties are consistent with the observer crossing the pulsar emission 
beam along the conal emission zone and the observer is likely to miss the core emission. Hence for PSR B1831— 04, the bright 
centrally located component can be interpreted as conal emission rather than core emission. Further there are situations where 
the PA below the core emission can deviate from the RVM model (see Mitra, Rankin & Gupta 2007), and hence the location 
of the infiexion point is difficult to determine. A careful examination using multi-frequency single pulse polarisation data 
is hence essential to find the location of the core emission. For our data set we have used the GL98 polarisation data sets 
wherever possible to find the location of the core emission which we discuss in the appendix. However incorrect identification 
of the core would clearly affect demonstration of the A/R effect. 

In several cases the central emission close to the core is weaker compared to the outer cones. Methods employed to find 
the peak of the core (like the Gaussian fitting scheme discussed in section |4]) have problems in establishing the location of the 
core component (like PSR B1738~08), and hence cause complication in discerning the A/R effect. 

Subpulse modulation (or drifting) provides another means to identify core and conal components in pulsars. The conal 
components tends to show subpulse drifting seen as a spectral feature in phase resolved fiuctuation spectra, while in the core 
emission this phenomenon is generally absent (Rankin 1986). It is sometimes possible that the central component of the pulsar 
shows drifting properties, in which case the central component cannot be core emission. We have verified our sample pulsars 
for such drifting of the central component based on the dataset of pulsar subpulse modulation published by Weltevrede et 
al. (2006, 2007). Three pulsars from our sample show central component modulation feature out of which two pulsars PSR 
B1508+55 and B1917+00 has a broad modulation feature which is consistent with core emission, while PSR B1944+17 shows 
clear drifting properties making this pulsar unsuitable for A/R analysis. 

3.3 Partial cone emission 

In several pulsars only a part of the emission cone is visible and these pulsars were identified as partial cones by Lyne & 
Manchester (1988). The PA inflexion point in these pulsars lies towards one edge of the profile. This kind of pulsar however 
poses a problem for finding the A/R effect, since estimation of the midpoint of the total intensity profile is incorrect due to 
missing emission. One such example in our sample is perhaps PSR B1859+03 as shown in Fig ([Sjl- The PA inflexion point 
coincides with the peak of the trailing component and has a sign changing circular, which is suggestive of core emission, 
however there is no emission component seen in the trailing part of the profile. Single pulse study by Mitra, Rankin & Sarala 
(2008) recently has shown that the missing partial cone regions occasionally flare up, and these flared proflles can be used to 
demonstrate the A/R effect. 
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Figure 1. B1857— 26 polarisation profile. The upper panel of the plot shows the intensity (solid line), linear polarisation (dotted line) 
and the circular polarisation (dashed line). The bottom panel shows the PA with error bars. This profile is an example where the core 
component, the inflection point and the sign changing circular are close to each other (see section 3 for details). 
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Figure 2. B1831— 04 polarisation profile (same as in Fig[TJ. Although this pulsar shows the sign change in circular it has a flat PA and 
a wide profile. Hence it is likely that the central bright emission component is not a core (see section 3 for details) 
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Figure 3. B1859-(-03 polarisation profile (same as in Fig[TJ. This pulsar has a clear sign changing circular indicating the location of the 
core emission. However the emission towards the trailing parts of the profile seems to be missing (see section 3 for details). 
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4 OBSERVATIONS AND DATA ANALYSIS 

For our analysis we selected profiles of 28 pulsars from three sources (see Table[l]): GMRT, AO and EPN. The GMRT (Swarup 
et al. 1991) operating at meter wavelengths is an array of 30 antennas with alt-az mounts, each of 45 meter diameter, spread 
out over a 25 km region located 80 km north of Pune, India. It is designed to operate at multiple frequencies (150, 235, 325, 
610 and 1000 — 1450 MHz respectively) having a maximum bandwidth of 32 MHz split into upper and lower sidebands of 16 
MHz each. It is primarily an aperture synthesis interferometer as well as a phased array instrument. The GMRT observations 
reported in this paper are done in the phased array mode of operation at RF 225 and 325 MHz, respectively, using the upper 
side band and a bandwidth of 16 MHz. In the phased array mode of the GMRT (Sirothia, S. 2000, Gupta et al. 2000), signals 
from any selected set of the antennas are added together after an appropriate delay and phase compensation to synthesise a 
single, larger antenna with a narrower beam. Signals of 16 MHz from each side band of each antenna are available across 256 
channels which can be summed in the GMRT array combiner. The summed voltage signal output from the phased array is 
then processed like any other single dish telescope signal, by the pulsar backend. At the frequencies of interest here, GMRT 
has two linearly polarised outputs for each antenna which are further mixed in a quadrature hybrid to produce two circular 
polarised signals which go further into the front end electronics and rest of the signal chain. The raw signals are finally post 
integrated by the pulsar backend and the data was recorded with the final sampling time of 0.512 msec. The time resolution 
for all the analysis done in this paper is the pulsar period divided by the sampling time. All the pulsars were observed for 
about 2000 pulse periods to obtain stable pulse profiles. The off-line analysis involved interference rejection, dedispersion and 
folding to obtain the average pulse profiles. 

The AO data set is from observations at the Arecibo Observatory, carried out between December 2001 and June 2002. 
The pulsars were observed at different frequencies using the following receivers : 430 MHz, L-band (1175 MHz), S-low (2250 
MHz), S-high (3500 MHz), C-band (4850-5000 MHz). The pulsar back-end used most frequently was the WAPP (Wideband 
Arecibo Pulsar Processor), though the PSPM (Penn State Pulsar Machine) was also used for some of the observations. The 
bandwidths of the observations were typically 25/50 MHz at 430 MHz and 100 MHz for all the higher frequency observations. 
The time resolution of the WAPP raw data was typically 128, 256 or 512 microsec. The data were analysed using the SIGPROC 
analysis package to get the average folded profiles. 

The data from the EPN archive were mostly obtained by GL98 using 76-m Lovell telescope at Jodrell Bank Observatory, 
except B1039— 19 at 4850 MHz, which was obtained by Kijak et al. (1997) using Effelsberg Observatory. Generally, resolution 
of EPN data is around 1°, which is much lower than that of GMRT or AO data. 

Our analysis followed these steps; (1) Identification of core components in profiles. We have used several established 
methods for discriminating the core components: (i) quasi-centrality of the component flanked by one or two pairs of conal 
components, (ii) associating the core component with the inflexion point of the polarisation position angle traverse and sense- 
reversing signature of circular polarisation under the core component (although the polarimetry was not always available 
in our data), (2) identification of conal pairs as leading and trailing components, primarily based on symmetric location of 
components on either side of the core one, (3) estimation of phase shifts according to eqn. ((5)1, and (4) estimation of the 
emission heights according to eqn. (|6]). This simple procedure is not always straightforward in practice. There are cases when 
components merge (especially as a function of the observing frequency) and their positions are not clearly determined. In 
these cases we decided to use a Gaussian fitting method introduced by Kramer et. al (1994). For consistency, we applied this 
method to all our profiles, even those in which all components were clearly distinguishable. The Gaussian fitting procedure 
was as follows: fix Gaussian parameters of identifiable components (position, amplitude, width) , put additional Gaussian with 
some initial parameters and allow the fitting method to find the best position, amplitudes and widths. For detailed information 
about Gaussian fitting method see Kramer et. al (1994). After applying this method to unresolved components we obtained 
a full and consistent set of component positions to be used in further analysis. 
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Table 1. List of pulsars chosen for analysis. The AO and GMRT observations are new, whereas the EPN profiles have been taken from 
the data archive http://www.mpifr-bonn.mpg.de/div/pulsar/data/ 



pulsar V (MHz) data source pulsar v (MHz) data source pulsar v (MHz) data source 



B0402+61 
B0410+69 
B0621-04 



B0656+14 
B0950+08 
B1039-19 



B1237+25 



B1508+55 



225 
408 
225 
408 
325 
410 
610 
1408 
445 
455 
225 
325 
1408 
4850 
157 
455 
1225 
2300 
3550 
4900 
5050 
225 
408 
610 
925 
1408 
1642 



GMRT 
EPN 
GMRT 
EPN 
GMRT 
EPN 
EPN 
EPN 
AO 
AO 
GMRT 
GMRT 
EPN 
EPN 
GMRT 
AO 
AO 
AO 
AO 
AO 
AO 
GMRT 
EPN 
EPN 
EPN 
EPN 
EPN 



B1541+09 



B1737+13 

B1738-08 
B1821+05 

B1831-04 



B1839+09 
B1848+08 
B1857-26 
B1859+03 



225 
325 
457 
1225 
2300 
4900 
225 
445 
457 
1225 
225 
325 
225 
445 
1225 
1460 
2200 
225 
325 
606 
2200 
445 
225 
1225 
2250 



GMRT 
GMRT 

AO 

AO 

AO 

AO 
GMRT 

AO 

AO 

AO 
GMRT 
GMRT 
GMRT 

AO 

AO 

AO 

AO 
GMRT 
GMRT 

EPN 

AO 

AO 
GMRT 

AO 

AO 



B1911+13 

B1916+14 

B1917+00 

B1944+17 

B1946+35 
B2002+31 

B2020+28 
B2053+21 
B2053+36 
B2315+21 

B2319+60 



325 
445 
1225 
445 
1225 
225 
325 
225 
325 
2300 
1225 
2300 
445 
445 
1225 
225 
1225 
225 
610 
925 
1408 
1642 



GMRT 

AO 

AO 

AO 

AO 
GMRT 
GMRT 
GMRT 
GMRT 

AO 

AO 

AO 

AO 

AO 

AO 
GMRT 

AO 
GMRT 

EPN 

EPN 
EPN 
EPN 



5 RESULTS 

Our preliminary data set was reduced to 23 pulsars, because it was hard to perform a reasonable fitting on 5 pulsars (B0950-I-08, 
B1859+03, B1946-I-35, B2020+28, B2315+21) due to either low S/N or problems with determining positions of core and/or 
conal components, even using a Gaussian fitting method. The determined positions of all components in a pulse profile, using 
the already mentioned Gaussian fitting procedure, were used for a further analysis. 

Firstly, we estimate the phase-shifts due to the A/R effect, and find the emission height tar using eqn. (|6]). For calculating 
the phase-shifts using eqn. ([Sjl, we first identify the leading and trailing component of a cone of emission using the location 
of the peak of the Gaussian fitted component as 0; and cf>t (this profile measure is often called as peak-to-peak component 
separation (PPCS), see for e.g. G&Ga and G&Gb, Dyks et al. 2004, Mitra & Li 2004). Correspondingly we denote the phase 
shift as A^ppcs- Wherever possible, /S.(j)ppcs was estimated for the inner and outer cones, separately. 

To estimate errors of single phase position measurements we used the formula introduced by Kijak & Gil (1997): 




where r is the observing resolution, oi is the root mean square of total intensity, and I is the measured signal level. Since W 
estimations need two measurements (left and right side of profile at particular intensity level) and Acj> estimations need three 
measurements ((^o, (j>\ and (j>t) we took 



(Tw = ^/ o-f + a? (8) 

'^A^ = ^J^l„+^l+^l (9) 

as an error for W and A(/!> respectively. 

The best example of a pulsar whose apparent components are very well reproduced by the Gaussian fitting method is 
PSR B1237-I-25 presented in Fig. (|4} at a number of frequencies. All B1237+25 profiles are fitted perfectly with 5 Gaussian 
components, being an ideal example of applicability of the Gaussian fitting technique used in this pulsar. This pulsar clearly 
shows the A/R effects at all frequencies, consistently over the spectrum from about 5 GHz to about 0.16 GHz (see also 
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B1237+25 




Figure 4. The profile evolution and fitted Gaussians of B1237+25 through a number of frequencies. This is a good example showing 
clear A/R effects at all frequencies. 



Table A. 5. 2 in Section A. 5 in the on-line Appendi^lj). The profiles are arranged in panels with ascending frequency. The 
original profiles are plotted with solid lines, while the fitted Gaussian components are drawn with dashed lines. The profiles 
are aligned with respect to the peak of the central Gaussian component identified as the core component. All the profiles of 
pulsar B1237+25 show a clear A/R shift. Srostlik & Rankin (2005) have found that PA steepest gradient traverse lags the 
sign changing circular polarisation zero point by about 0.4°. However if one takes this value into account, B1237+25 still 
shows the A/R shift. 

Our main results, which will be discussed later on, show the commonness of A/R effects in pulsar profiles. Within our 
sample of 23 pulsars 7 show clear A/R effects, 12 of them (doubtful cases) show a clear tendency towards this effect, while 
the remaining 4 are so-called counter examples. When it comes to single A<j!> measurements we found that 45 measurements 
show negative shifts, which is expected for A/R effects, 23 measurements are doubtful within uncertainties, while remaining 
19 show opposite shifts. However taking into account only the outer cones (where more than one cone is seen) reveals that 
30 shift measurements show A/R effects, 18 are doubtful cases whereas 14 show no A/R effects at all. For the majority of 
non-A/R cases, the problem lies mostly in misidentifying core components. Some of the cases are partial cones which also 
makes it difficult to measure A/R shift properly. For a detailed discussion of various issues regarding A/R calculations see 
section ([3|. For a discussion of the results see section (|6}. 

Note that all individual results are available in the form of figures and detailed tables in the electronic form as an appendix 
to this paper. 



^ http: / /astro. ia.uz.zgora.pl/~chriss/Krzeszowski_2009_appendix.ps 
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6 CONCLUSIONS AND DISCUSSION 

Observationally, the peak of the core component lagging the centre of the outer conal component has always been considered 
as an evidence for A/R effects operating in pulsars. Those effects have been demonstrated in a few pulsars so far: PSR 
B0329+54 by Malov & Suleimanova (1998) and G&Ga between 103 MHz to 10.5 GHz and PSR's B0450-18, B1237+25, 
B1821+05, B1857-26, B2045-16 and PSR B2111+46 at 325 MHz by G&Gb. In this work, our aim was to verify the A/R 
effects for a large sample of pulsars. While our sample consists of 23 pulsars, we note that our analysis is aided by using 
multi-frequency measurements, where each frequency serves as an independent data set. Hence, although our sample consists 
of three pulsars common with G&Gb, the measurements are done at different frequencies and can be treated as independent 
check for their results. In our data set of 62 pulse profiles measured, the A/R effects are clearly seen in 30 profiles. Fig. (O 
shows the distribution of the measured phase shifts Aifippcs for our sample of pulsars (for the data see the Appendix). The 
tendency for the negative shift, which is the signature for A/R effect, is clearly visible. 

Out of the four counter-examples noted in our sample, the most outstanding one is PSR B1831— 04 (Fig. [2|. We think 
that the central component in this pulsar might be incorrectly identified as the core emission (as discussed in section 3). The 
flat PA traverse seems to suggest that the line of sight is not cutting the overall beam centrally (see another example of broad 
profile pulsar B0826— 34 discussed by Gupta, Gil, Kijak et al. 2004), so the core emission might be missing. If this is the case 
then the central component could arise from grazing the third, innermost cone of emission, although its non-central location 
is a concern in this model. However, if our hypothetical 3-rd cone interpretation is incorrect, then B1831— 04 remains the 
strongest counter-example of A/R effects. 

The three other counter-examples B1738— 08, B1848+04 and B1859+03 show a marginal deviation from the A/R effects 
and are probably also related to our inability of identifying the location of the core component. In the first case (B1738— 08) 
the core component is not apparent, although a blind Gaussian fitting method indicates a central "core" component (see Fig. 
A. 9 in the on-line Appendix). The second pulsar (B1848-I-04) is again a typical example of a broad profile pulsar (see Fig. 
A. 13 in the on-line Appendix), where the core component is likely to be missed. The third pulsar (B1859-I-03) (see Fig. A. 4 
in the on-line Appendix) is a clear example of a pulsars with missing component (in this case the trailing one). Therefore, 
we conclude, that all 4 cases of strong counter-examples reported are the result of mis-identification of the core component. 
Sometimes weak conal emission or merged conal components leads to wrong identification of the emission components. For 
example in PSR B1541+09 the leading conal component at 1225 and 2330 MHz is difhcult to model. In cases like PSR 
B1944+17 and B1737+13 the weak conal emission poses difficulties in determining the location of the components. It is likely 
that the same set of problems occur in the group of doubtful cases (i.e. pulsars which show A/R effects at some frequencies 
but not in others). This strengthens our general conclusions that A/R effects are apparent in pulsar profiles. 

In cases where A/R effects are seen between the core and conal components in pulsars, implies that the core emission 
necessarily has to originate closer to the stellar surface than the conal emission (see for e.g. G&Ga, Dyks et al. 2004). The delay 
emission heights obtained by eqn. ^ are essentially the conal emission heights with respect to the core emission. Alternatively 
geometrical methods can be used to find the conal emission height tkg as specified by eqn. Q. A comparison between the 
conal emission heights measured by these two methods is shown in Fig. ([6]). This plot can be compared to Fig. (29) of BCW 
where a similar effect is noted, although the delay heights there were measured from polarisation observations and they used 
log-log axes. In this figure the clear A/R cases are marked with filled circles, doubtful cases are marked with open circles and 
the so called counter-examples are marked with stars. Note that many of the negative tar values, in Fig. (|6]) correspond to 
the counter-examples in our sample. It is noteworthy that statistically the relation vkg > '"ah, seems to hold well. There are 
primarily two reasons which can explain this effect. Firstly the tar values can be underestimated if the putative core emission 
arises from a definite height above the polar cap. In order to square the two height estimates by making the cluster of points 
in Fig. ((61 to be located on the 45° line, a nominal average core emission height of about 100 km is needed. Alternatively, the 
geometrical emission height rxG might be overestimated, since the assumption prevailing in estimating these heights is that 
the measured widths entail the last open field lines, i.e. s — 1 (see eqn. [2)). If the parameter s was around 0.7, the two height 
estimates would square with each other. As discussed by Mitra & Li (2004), a careful and systematic study of the relation of 
the core emission w.r.t. the polarisation position angle traverse might help to resolve this issue. 
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Figure 5. Distribution of measured phase shifts Acj>ppcg. Note that a clear tendency for negative shifts is seen in the sample. 




Figure 6. Relation between the delay emission height tar given by eqn. (|6]l and the semi-empirical geometrical height estimate rxG 
given by eqn. l(4)l. Note that statistically rxG > '"AR as seen by the excess of points lying above the tar = f'KG li^G in the figure (see 
legend to distinguish different groups of pulsars). See text for details. 
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